Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  933 1559 (N39W11 L=163)   676    178        1 ( 879/1485- 999/1615)  17.1 W05-W18 N34-N41
      2 1213 1536 (N28W28 L=179)   351     53        0 (1192/1507-1236/1565)  76.0 W26-W29 N26-N30
      3 1349 1510 (N23W34 L=186)   962    175        3 (1267/1470-1430/1539) -37.9 W29-W38 N18-N25
      4  293 1208 (N34E45 L=106)   313     85        2 ( 264/1188- 323/1229) -18.4 E48-E42 N32-N35
      5  463 1156 (N27E32 L=120)   517    153        1 ( 402/1135- 529/1182) -11.3 E36-E26 N25-N29
      6  691 1086 (N17E18 L=134)   608    119        1 ( 632/1067- 758/1115) -35.5 E20-E13 N15-N20
      7 1335 1058 (S01W18 L=169)   978    211        2 (1229/1025-1440/1101) -26.6 W11-W23 S05-N01
      8 1906 1048 (S20W61 L=213)   510     22        0 (1889/1027-1923/1069) -90.0 W59-W63 S21-S18
      9 1300  981 (S04W14 L=165)   264     98        0 (1255/ 953-1352/ 996)  -8.1 W10-W17 S05-S03
     10 1517  896 (S15W25 L=176)   504    100        1 (1472/ 867-1564/ 931)   0.0 W21-W28 S16-S13
     11 1232  734 (S16W04 L=155)   760    193        1 (1199/ 659-1276/ 840)  79.7 W02-W06 S21-S10
     12  912  593 (S16E18 L=133)   649    150        1 ( 816/ 561- 974/ 623) -48.8 E22-E15 S19-S11
     13  629  440 (S18E42 L=109)   269     61        1 ( 602/ 418- 647/ 465)  36.9 E43-E39 S19-S16
     14  598  392 (S20E46 L=105)  1050    258        5 ( 514/ 305- 668/ 463) -81.3 E49-E44 S27-S14
     15  521  383 (S19E53 L= 98)   376     72        2 ( 498/ 357- 547/ 406)  63.4 E55-E50 S21-S17
     16  902  168 (S46E43 L=108)   863    138        1 ( 883/ 129- 923/ 208)  41.6 E48-E39 S49-S41
  

    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