Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  815 1616 (N22E29 L=240)   579    131        0 ( 757/1589- 882/1640)  45.0 E31-E25 N19-N25
      2  353 1582 (N11E63 L=206)   585     66        1 ( 329/1559- 389/1603)  36.9 E66-E59 N09-N12
      3 1676 1428 (N35W33 L=303)  1347    263        0 (1562/1354-1766/1493)  -6.0 W22-W43 N33-N36
      4  335 1373 (S02E53 L=216)   849    186        2 ( 282/1312- 405/1433)  49.4 E57-E49 S05-N02
      5  603 1274 (S03E30 L=239)   220     35        1 ( 584/1247- 623/1295)  90.0 E31-E28 S04-S01
      6  326 1234 (S10E50 L=219)   538     87        2 ( 288/1209- 363/1265) -15.9 E53-E46 S11-S08
      7  927 1191 (S00E09 L=260)   198     82        0 ( 911/1151- 947/1232) -36.9 E11-E07 S02-N01
      8  556 1164 (S10E31 L=238)   432    121        0 ( 538/1103- 592/1216) -45.0 E32-E27 S12-S06
      9 1027 1129 (S02E02 L=267)   207     56        1 (1006/1108-1057/1146)  45.0 E03-E01 S02-N00
     10 1916 1096 (N23W63 L=332)   720     98        0 (1901/1046-1930/1144)  80.5 W60-W65 N19-N25
     11 1848 1065 (N18W54 L=323)   360     39        1 (1826/1046-1868/1092)  71.6 W51-W56 N16-N19
     12 1174  684 (S22W18 L=287)  1108    349        4 (1132/ 587-1244/ 833) -71.6 W13-W22 S27-S12
     13  702  501 (S46E07 L=262)  1594    402        3 ( 553/ 432- 880/ 573)   3.9 E21-W08 S48-S43
  

    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