Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  990 1727 (N46E18 L=145)   329     48        0 ( 958/1714-1018/1745)  11.3 E20-E15 N44-N47
      2 1138 1648 (N42E04 L=160)   417    104        0 (1084/1626-1181/1677) -14.0 E08-W00 N41-N43
      3 1063 1597 (N37E08 L=155)   242     48        1 (1041/1576-1093/1615)  26.6 E09-E05 N35-N38
      4  339 1447 (N15E57 L=106)  2730    418        8 ( 273/1298- 399/1601) -84.6 E62-E53 N04-N25
      5  532 1313 (N10E36 L=127)   525    126        2 ( 496/1266- 572/1350) -15.9 E39-E32 N06-N11
      6 1063 1264 (N15E01 L=162)   259     63        0 (1036/1232-1086/1290) -33.7 E03-W00 N13-N16
      7  644 1192 (N04E26 L=137)   266     81        0 ( 627/1153- 665/1236) -59.0 E27-E24 N01-N06
      8  562  964 (S11E28 L=135)  1046    246        3 ( 532/ 857- 612/1096) -70.3 E30-E23 S16-S02
      9 1254  964 (S00W15 L=178)   457    118        0 (1196/ 931-1318/1001)  29.7 W11-W18 S02-N02
     10 1778  980 (N09W53 L=216)   311     51        1 (1758/ 955-1796/1000)  53.1 W51-W54 N06-N10
     11 1191  843 (S08W13 L=176)   221     81        1 (1167/ 802-1213/ 888)  80.5 W11-W13 S11-S05
     12  390  718 (S29E40 L=123)   280     67        0 ( 369/ 694- 415/ 744) -26.6 E42-E38 S30-S27
     13 1464  665 (S15W35 L=198)   308    129        0 (1425/ 619-1519/ 689)  45.0 W32-W38 S18-S13
     14  225  639 (S38E59 L=104)   690    181        2 ( 203/ 610- 244/ 677)  50.2 E62-E55 S40-S34
     15 1330  552 (S24W29 L=192)   259     63        0 (1296/ 537-1366/ 568)  26.6 W26-W30 S25-S22
     16 1271  421 (S34W30 L=193)   306     68        0 (1241/ 401-1308/ 446)  -8.1 W26-W33 S35-S32
     17 1540  377 (S32W53 L=216)  1442    285        1 (1462/ 336-1659/ 432)  29.7 W45-W61 S35-S26
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y